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Abstract 

We calculate o?o? s order QED contribution to J/tp production in the pp(p) — ► J/tp + cc color- 
singlet process at the Tevatron and LHC in the framework of nonrelativistic QCD. The contribution 
of the interference between the a 2 a 2 s QED and QCD is also taken into account. The J/tp 
production associated with a charm-quark pair could be a measurable signal at hadron collider 
experiment. Our calculations show that by including the QED contribution, the pt distribution is 
enhanced by a factor of 1.5 (1.9) at the Tevatron (LHC) at pt = 50(100) GeV. In addition, the 
polarization of J/tp turns from unpolarized in all region to increasingly transverse when pt becomes 
larger. 

PACS numbers: 12.38.Bx, 14.40.Gx, 13.85.Ni 
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I. INTRODUCTION 



Heavy-quarkonium production provides important tests on our understanding of QCD 
in both perturbative and nonperturbative aspects. Since the heavy quark mass is larger 
than Aqcd, heavy quarkonia are commonly thought to be nonrelativistic bound states. 
Now the nonrelativistic QCD (NRQCD) effective theory[l] is widely accepted to study the 
heavy-quarkonium phenomenology (for a review see^). In the framework of NRQCD, the 
production rates of heavy quarkonium could be calculated with a rigorous factorization 
formula based on the double expansions of the QCD coupling constant a s and heavy quark 
relative velocity v in the heavy mesons. 

The vital difference between NRQCD and traditional color-singlet model is that the 
NRQCD allows the production of heavy quark pair in color-octet configuration over short 
distance, and the color-octet state consequently evolves into the physical heavy quarkonium 
through the nonperturbative emission of soft gluons. The introduction of the color-octet 
mechanism has successfully reconciled the orders of magnitude discrepancies between 
the experimental data measured by the CDF Collaboration 4J and the color-singlet model 
theoretical predictions Q]. However, the CDF former 6] and latest measurements of J hp 
being almost unpolarized conflict with the NRQCD transverse polarization predictions |8|. 
So the J/if> production mechanism is still a big challenge, and various theoretical attempts 



can be found in Refs 



111- 



Recently, substantial progress has been achieved in the calculation of high order QCD 
corrections to J /if) hadroproduction, which are very helpful to clarify the big puzzle. The 
next-to-leading-order (NLO) QCD corrections [12J to color-singlet gg — > J/ipg production 
enhance the leading-order (LO) cross section by a factor of 2. Moreover, the NLO result has 
a large impact on the transverse momentum (p t ) distribution of the production rate. It is 
because the p t distribution of the LO result behaves as 1/pj? , while the NLO result behaves 
as 1/pf due to involving the new topological Feynman diagrams. It is also found Jij, 14] that 
the color-singlet J / if) production in association with the cc pair process gg — > J /if) + cc also 
has sizable contribution to the a\ (NLO) corrections. This process is also investigated in the 
kx factorization formula [l5|. Furthermore, it is reported 16] that the polarization of J /if) 
production via the color-singlet 3 Si state turns from more transverse to more longitudinal 
after including the NLO QCD corrections, and most of the J /if) produced via color-octet 1 Sq 
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and 3 Si states at QCD NLO are still almost in a transverse polarization state [it| in the large 
Pt region. Toward the resolution of the puzzle on J/ip production mechanism, the study 
has been extended to looking for other possible experimental observable J/ip production 
processes. J/ip production associated with a photon in the process pp(p) — > J/ip + 7 is 



studied in Ref. [18|. Its NLO correction has recently been performed in Ref. 

ysis has been obtained in Ref. 



2l| 



19ll and the 
3- Double 



real next-to-next-to-leading-order (NNLO) ana 
J/ip hadroproduction has been studied in Ref. 

Among these attempts, the study on the J/ip associated production, i.e. J/ip + cc + X 
is very interesting. The reasons come from two sides. Experimentally, the cross section 
of the associated production could be measured by searching for J/ip together with DD 
pair, where D denotes all possible charm hadrons such as D , D + , D*, A c , etc. Theoretically, 
it is noticed[22] that at NNLO the infrared divergence appears due to the soft interaction 
between the associated core quark and the cc pair in charmonium. Hence such processes are 
nonfactorizable at NNLO. And the color-transfer enhancement is introduced in the double- 
heavy-quark-pair production process. One persuasive example is the recent measurements in 
e + e~ annihilation at center-of-mass energy yfs = 10.6 GeV at B factories for the inclusive 
and exclusive charmonium production via double charm-quark pairs. Belle collaboration 
first found 123 ] 

rr{< + < ~ — J/< ' - re) ... , ., , , 



J/ip + cc) 



0.59;°;^ ±0.12, 



241 ] . The experimental result is about 6 times 



cr(e + e — > J/ip + 
which is confirmed by their latest analysis 

larger than the NRQCD predictions at LO in a s and v [25|]. The problem is resolved by 

the recent NLO QCD corrections to inclusive J/ ip production at B factories 26] . Another 

example is that the associated production process J/ip + cc is also found to be important in 

J/ip photoproduction [271] . 

During the procedure of studying the charmonium production in e + e~ annihilation, it is 

found that the QED contribution is sizable since the J PC = \— state J/ip can be produced 

via one photon fragmentation. For example, __(a) the cross section is very large in J j ip 

(b) in e + e~ — > J/ip + r\ c process, the 



electromagnetic production at e + e~ colliders 
contribution of the interference between QCD and QED could enhance that of the QCD 
result by 20%, although the ratios of the QED contributionjxi the QCD contribution is 
about 1% 



n 



29j; (c) the production rate of e + e — ► J/ip + J/-?/>|30! at LO is about 3.7 times 
larger than that of e + e~ — » J/ip + 7/ c [29j at t/s = 10.6 GeV, although the QCD correction 



to the former process is negative 3l| ; (d) when the center-of-mass energy ^/s > 20 GeV 32] , 
the J/if> + cc produced in e + e~ annihilation through two virtual photons will prevail over 
that through one virtual photon. 

In this situation, it is necessary to examine what will happen in J/ip hadropro duct ion 
after including the QED contribution. Our previous work 33J has shown that the photon 
fragmentation process up to QCD NLO is very important for the p t distribution of the 
production and polarization for J/if) hadroproduction, though its contribution to the total 
cross section is negligible compared to the QCD contribution. As mentioned above, the J/ip 
associated production is also a very interesting process. Following our previous work, we 
study in detail the QED effect on J/if> + cc hadroproduction at the Tevatron and LHC in 
this paper. The rest of the paper is organized as follows: In Sec. II, we briefly introduce 
the basic formulas and input parameters used in our calculation. In Sec. Ill, we discuss 
the QED contribution to the p t distributions of J/if> production and polarizations for the 
associated J /if} hadroproduction at the Tevatron and LHC. The conclusion is presented in 
Sec. IV. 



II. BASIC FORMULAS 

According to the NRQCD factorization approach[l|], the J /if) + cc production rate in 
hadron-hadron collisions could be expressed as 

a\pp ->J/V + c + c + X] = ^P j dx 1 dx 2 G i i v Gji v x fr[i + j -> (cc) n + c + c + X] (0%), (2) 



where p is either a proton or antiproton. The short distance part a represents the partonic 
production of cc with quantum number n. (O^) is the nonperturbative long-distance matrix 
element that parametrizes the transition of the cc pair into J/if>. 

The QCD contribution to pp(p) — > J/ip + cc + X starts at a* order with two independent 
partonic processes: gg — > cc[ 2S+1 Lj, 1(8)] + c + c) and qq cc[ 2S+1 Lj, 1(8)] +c + c, where q 
represents all possible light quarks u,d,s. cc[ 2S+1 Lj, 1(8)] denotes the specific state of cc with 
total spin S, orbital angular momentum L, total angular momentum J and 1 (color-singlet) 
or 8 (color-octet). The contributions of J/ip production via the color-singlet cc[ 3 Si,l], 
color-octet c c\ 3 S \, 8], and ccf 1 ^,^] states in the gg fusion processes have been considered in 



Refs 



.0 
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34]- Since only the cc[ Si, 1] state can directly couple with one photon, which 



may bring kinematic enhancement, in this work, we consider the QED contributions from 
two processes: 

g(ki)g(k 2 ) -> cc[ 3 S u lj(p 1 ) + c(p 2 ) + c(ps), 

^1)^2) - cc[ 3 S 1; l]( Pl ) + c(p 2 ) + c(p 3 ). (3) 
The amplitude at parton level is expressed as 

M [i(h)j(k 2 ) - cc[ 3 S 1; U( Pl ) + c(p 2 ) + c(ps)] = £ X>i5 s a|l^) 

si,S2 3fc,3Z 

(3*; 3Z|1) x [i(A:i)j(fe) -> c fc (^, si) + c,( j, s 2 ) + c(p 2 ) + c(p 3 )] , (4) 

where z,j are gluon or quarks, (3/c;3Z|l) = 5m/\/N^ and (si; s 2 |lS' z ) are S77(3)-color and 
.St/ (2)-spin Clebsch-Gordon coefficients for cc pair projecting onto the color-singlet 3 Si state. 
At LO in v, the projection operator of Dirac spinor is 

p i,sM = J2( Sl] s 2|l^)^(y, s 2 ) = ——ftS^fa + 2m c ). (5) 

Sl,S 2 _ _ V 

To obtain the polarization distribution, the polarization vectors e M (A) are kept explicitly 
during our calculation. The partonic differential cross section for polarized J/ip can be 
written as 

^ = Ae(X) ■ e*(A) + £ A ijPi • e(A) B ■ e*(A), (6) 

where A = Tx,T 2 ,L. e(Ti), e(T 2 ) and e(L) are the two transverse polarization vectors and 
the longitudinal one for J/ip, respectively. A and Ajj are the coefficients. The polarization 
parameter a is defined as 

/ \ = da T /dp t -2da L /dp t 
a[Pt) da T /dp t + 2da L /d Pt - U 
We chose physical gauge for gluons to avoid computing the ghost diagrams. The gauge 

invariance is checked to ensure the validity of our results by replacing the gluon polarization 

vector with its momentum in numerical computation. Moreover, the long-distance matrix 

element (O(^) is estimated from the leptonic decay of J/ip using the relation 

^ = 2^(27+1)1^(0)1^ 
4ir 

At NLO in a s and LO in v, we have T(J/i/j -> e+e_] = ^(1 - 16a s /37r) |i2^(0) | 2 . We 
employ the Feynman Diagram Calculation package 36( to generate the Feynman diagrams 
and calculate the invariant amplitudes numerically with the basic formulas mentioned above. 



The pp(p) — > J/ip + cc process is part of the QCD real corrections to the a s a 2 order J/if> 
inclusive production. So we chose the same set of numerical inputs as our last work with 

• QCD coupling constant: a s (M z ) = 0.118; 

• renormalization scale /i r and factorization scale [if. \i r — fif = fi = 

• QED coupling constant: a = 1/128; 

• charm-quark mass: m c = 1.5 GeV; 

• Jfij) leptonic decay width: r(J/ip — > e + e~~) = 5.55keV 371 ]; 

• long-distance matrix element (0"[ ) = 1.35GeV 3 from the leptonic decay width; 



PDF set: CTEQ6M|38j; 



• kinematic cut: p t > 3GeV; and \yj/^\ < 0.6 (Tevatron), \yj/$\ < 3.0 (LHC). 

III. QED CONTRIBUTION TO pp(p) -» J/V> + cc 

To investigate the QED effect in pp(p) — > J/ip + cc, we present the calculation in three 
parts : QCD, QED and their interference. The QCD part of gg — ► cc[ 3 Si, 1] + c + c 
includes 42 Feynman diagrams shown in Fig. [lj While there are only 38 Feynman diagrams 
for the a 2 a 2 s QED part, which are divided into two types and each type forms a gauge 
invariant subset. The Type I diagrams shown in Fig. [2] can be obtained by replacing the 
gluon lines between quarks in the QCD diagrams by photon lines. The Type II diagrams 
shown in Fig. [3] are similar with some of the gluon fragmentation diagrams for color-octet 
cc[ 3 Si,8] + c + c production when 7* — > cc[ 3 Si, 1] is changed into g* — > cc[ 3 Si,8]. So 
the partonic differential cross sections of the two types QED diagrams have the same 1/pf 



scale behaviors as the QCD ones 341] . Then the contribution of Type I diagrams will be 
suppressed by (a/a s ) 2 compared to the color-singlet QCD result in both the small and large 
p t region. However, there is a kinematic enhancement in the Type II diagrams. The reason 
lies simplify in the fact that the virtuality of the photon in Type II diagrams are fixed to 
4m 2 ., whereas it varies from Am 2 to p 2 order in Type I diagrams. When p t is large enough, 
the kinematic enhancement may compensate the suppression factor (a/a s ) 2 . Furthermore, 
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the value of a s (//o) will become smaller as pt increasing. We verified that the contribution 
of Type I diagrams is dominant in small pt region and drops fast when p t increases. For the 
qq process, there are 7 QCD and 16 QED Feynman diagrams shown in Fig. m and the QED 
Feynman diagrams also can be divided into two types. 




FIG. 1: The typical Feynman diagrams for the QCD part of gg — > cc[ 3 Si, 1] + c + c. 
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FIG. 3: The typical Type II Feynman diagrams of photon fragmentation in the LO {a 2 a 2 s ) QED 
part of gg — ► cc[ 3 Si, 1] + c + c. 











FIG. 4: The typical Feynman diagrams of the LO process of — > cc[ 3 Si, 1] + c + c for QCD and 
QED. 



For the QCD contribution alone, we calculate the total cross sections and pt distributions 
of J /if) production (dashed lines in Fig. [5] (a)) and polarizations (the dashed lines in Fig. [6) 
for the Tevatron and LHC The results are in agreement with the early predictions 13l.ll4j.l34] 
when the difference of the parameters are taken into account. And the corresponding total 
cross sections are 

11 qcd ipP ~^ J /ip + cc + X) = 8.7 nb (Tevatron), (9a) 
a QCD (pp-^ J/tp + cc + X) = 169.6 nb (LHC). (9b) 
For the QED contribution alone, the total cross sections are 



17 qed (PP ~^ J/ip + cc + X) = 52.3 pb (Tevatron), 
°qed(pP^ J/^ + cc + X) = 1.13 nb (LHC). 



(10a) 
(10b) 
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FIG. 5: (a): The pt distributions of the production rate for J/tp + c + c hadroproduction 
at the Tevatron (^/i =1.96TeV, \y\ <0.6, 3GeV<pt<50GeV) and LHC (y/s =14TeV, \y\ <3, 
3GeV<pt<100GeV). (b): the pt distributions of the "total" with /j, r = [if = pL dependence. Where 
the total refers to the sum of the QCD, QED and their interference results. 

Since the integrated cross sections are mostly coming from the contribution in small pt 
region, the total cross sections of the QED part are about (a/a s ) 2 ~ 1/100 times smaller 
than the QCD ones for both the Tevatron and LHC. 

The p t distribution of J/ip QED production is shown in Fig. [5] with lower (upper) dotted 
line for the Tevatorn (LHC). It can be seen that the QED contribution is about 2 orders 
of magnitude less than the QCD contribution at pt — 3 GeV and the ratio of the QED 
contribution to the QCD one is 0.4 (0.7) at p t = 50 (100) GeV for the Tevatron (LHC). It 
means that the QED contribution is indeed negligible in the small p t region, but is important 
in the large p t region, especially at the LHC. The p t distributions of J/ip polarization from 
the QED part for the Tevatron and LHC are plotted with the dotted-dashed lines in Fig. [61 
Both of the curves show that J/ip is unpolarized in the small pt region and becomes increas- 
ingly transverse in the region of large pt- It is because the QED part is dominated by the 
contribution from Type I diagrams in the small pt region, while the Type II fragmentation 
contribution dominates in the large p t region. 

For the QED part alone, the contribution is comparable to the QCD one in the large p t 
region and particularly the polarization behavior of the former is very different from that of 
the latter. Therefore, we also calculate the interference between QED and QCD parts. It is 
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FIG. 6: The pt distribution of polarization parameter a for the J/ip + c + c hadroproduction at 
the Tevatron (a) and LHC (b). Here, the total refers to the sum of the QCD, QED and their 
interference results. 

found that the interference between the QED and QCD parts of pp(p) — > cc[ 3 Si, 1] + cc is 
positive and the relative phase angle is close to n/2. The total results are shown in Fig. [5] 
(b). The theoretical uncertainties of the p t distribution for J/ip production are studied 
by changing the renormalization and factorization scale fif = fi r from /xo/2 to 2/j,q. They 
correspond to the uncertainty bands of J/ip pt distribution shown in Fig. [5] (b). The pt 
distribution of J/ip production is enhanced by 50% compared to the QCD result at p t = 50 
GeV at the Tevatron. At the LHC, the enhancement is 190% at p t = 100 GeV. The 
Pt distributions of J / ip polarization are shown in Fig. El It is shown that the polarization 
parameter a(p t ) for the QCD part is almost equals to zero in all regions at both the Tevatron 
and LHC. But it becomes increasing slowly (see the solid line in Fig. [6]) when QED part is 
included in the full calculation. 



IV. CONCLUSION 

Until now, people still have not found a convincing mechanism to explain the J/ip 
hadroproduction at the Tevatron. It is suggested [3] that the associated production channel 
J/ip + cc, which may be observed through measurement of the J/ip together with at least 
one charmed hadron, may provide further insight to the mechanism responsible for the J / ip 
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hadroproduction. Both the color-singlet and color-octet contributions for the associated J/ if) 



production have been calculated at at A s order in Ref. 34j. In this work, we have presented the 
a 2 a 2 s order calculation for the associated J /if) production at the Tevatron and LHC Our 
calculation also include the interference terms between the a 2 a 2 s order QED part and a A 
QCD color-singlet part. Our results indicate the QED photon fragmentation effect is very 
important, and it has a large impact on p t distribution of J/if> production and polarization 
in large pt region, especially for the LHC. 

The QED part for pp(p) - W + cc is also part of the NLO QCD corrections to photon 
fragmentation inclusive J / if) production considered in our previous work [33( . Both of them 
scale like so unlike the J/if> QCD production case, the impact of the QED J /if) associ- 
ated production contribution on the photon fragmentation inclusive J /if) production is not 
very large for the pt distribution of J/if> production and polarization. 

In this work, we do not study the T production associated with the bb pair, though T 
can also be produced by photon fragmentation. There are two reasons. One is the electric 
charge of bottom-quark is only half of that of charm-quark. It will provide an additional 
(1/2) 4 suppression factor compared to the charm-quark case. The other is the mass of T, 
which is about 3 times larger than the J /if) mass, which makes the photon fragmentation 
effect for T not as outstanding as that for J/if> case. 

This work is supported by the National Natural Science Foundation of China under 
Contract No. 10775141 and by the Chinese Academy of Science under Project No. KJCX3- 
SYW-N2. 
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